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1.0 INTRODUCTION 


The Airborne Laser Ranging System has been proposed as an effort 
whereby measurements of the earth's crustal shift can be detected to 
within centimeters over a region of approximately 200 x 400 km. Hope- 
fully by gaining an understanding of this shift, scientists will be able 
to more accurately comprehend and predict earthquakes. 

Previous analyses have employed spaceborne laser ranging systems. 
Because of obvious problems arising in the use of spacecraft, this task 
is devoted to determining the requirements and limitations of employing 
an airborne system. 

Essentially the system consists of an aircraft which flies over a 
grid of ground deployed retroreflectors. The aircraft flies at two alti- 
tudes and presently makes six passes over the grid. 

The retroreflector baseline errors are assumed to result from 
measurement noise, a priori errors on the aircraft and retroreflector 
positions, tropospheric refraction, and sensor biases. 

1*1 Experiment Description 

Figure 1.1-1 shows the configuration of the present grid system of 
retroreflectors. There are currently 15 retros in an area of 40 x 80 km. 

The aircraft position is determined by interpolating between the 
position of the aircraft at the beginning of a pass over the grid and the 
end of the pass. The aircraft passes start and end 5 km on either side 
of the grid and the A/C usually flies in the lengthwise direction of the 
grid at a constant speed of 800 km/hr. Pulses are attempted at a con- 
stant time interval. 

At each pulse time, a test is performed to select all of the sta- 
tions in view (those above the 20° elevation cutoff angle), and of 
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these, the stations upon which measurements will be taken, are chosen. 

At this point, the sensitivity matrix is computed and the Kalman filter 
updates are performed. 

Thi:§ process repeats until all passes have been completed. Then the 
RSS errors are computed in cartesian and baseline coordinates between 
each station and the master station, which is currently station 1. If at 
any time, not enough stations are visible, the aircraft simply moves on 
to the next time. 


,11.12.13.X4.1S 
,7 ,e ,9 .10 
1.2 3 4 5 

20 km . 


Figure 1.1-1 
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2.0 Error Analysis 

2.1 Error Models 


2.1.1 A Priori Station Errors and Constraints 


The station positions are assumed to be known in an earth-fixed 
coordinate system to a level of 100 m (l^J), independent errors in each 
coordinate direction. 

Because the multiple laser measurements from aircraft (geometric 
measurements) are incapable of determining the global location of this 
grid, certain station coordinates are assumed to be perfectly known, thus 
providing the basis for a well-posed problem. The station 1 position is 
assumed, thus holding the grid in translation. Then the y and z 
coordinates of station 5 are fixed. This fixes the direction (though not 
the length) of the 1-5 baseline. In effect, the x-axis of the coordinate 
system has been defined. Finally, the z component of station 13 is 
assumed known. These mathematical constraints do not affect the inter- 
station baseline lengths. 

2.1.2 A Priori Aircraft Errors 


It is assumed that the A/C position is a gaussian random variable 
with a standard deviation of 6.3 m in each coordinate .direction. 

Currently, all runs are made with this error holding every time a 
pulse is taken. A procedure has been implemented in the Kalman filter 
for estimating A/C position and velocity in the presence of process noise 
(see Appendix A). This option, which could significantly improve base- 
line estimation accuracies at realistic sampling intervals (^^0.2 sec) has 
not been exercised because it has not been possible to obtain realistic 
estimates of the process noise variance and time constants. 
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2.1.3 Measurement Noise and Bias 

The laser range measurements are assumed to be corrupted by a bias 
and a mean-zero, white sequence. The bias is selected from a mean-zero 
distribution having standard deviation of 1 cm. Each beam has an 
independent bias and these are solved for in the Kalman filter. The 
white sequences added to the range measurements are independent on each 
beam and are gaussian with standard deviation o^ 1 cm. 

2.1.4 Refraction Errors 


The model for tropospheric refraction is summarized in [Error 
Analysis of the Spacelab Geodynamics Laser Ranging System, B. P. Gibbs 
and E. M. Haley, BTS-TR-78-52, February, 1978]. It appears that the 
dominant cause of residual refraction error (after range correction for 
refraction) is errors 1n the pressure and the gradient of PTK at the 
retroref lector sites. The procedure hypothesized is to take readings of 
temperature and pressure at various sites in and around the grid and then 
fit this data with a polynomial fit to pressure and PTK. This has been 
done with standard deviations of 1.4® on temperature measurements and 1 
mbar on pressure measurements at each of 15 meteorological sites. This 
data was used to fit pressure and PTK with quadratic polynomials (6 
coefficients in each approximation) modeling the horizontal variation and 
a linear polynomial modeling the horizontal variation of the altitude 
dependence. 

This regression leads to a covariance matrix for errors at the 
retroref lectors and this is used in the error analysis. 
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2.2 Se quential Error Analysis 

Define; 

6a » x(0|0) - x(0) , total error in the prior estimate. 

6x(k|ji) = x(k|)^) - x(k) k > A . 

6p = p - p , total error in a set of parameters affecting 
the measurements. 

fiX^Ckja) = S^(k|A)6a: Error in estimate caused by error in 

prior estimate. 

6Xy(k|A) = Sy(klA)v: Error in estimate caused by measurement 

noise. 

“ Sw(kU)w: Error in estimate caused by process noise. 

* 

6Xp(kjt) = Sp(kl.t)5p; .Error in estimate caused by parameter 
errors. 

fix(klt) = 6Xg + 6Xy + 6X^ + 6Xp 

x(k|t): Error in estimate caused by modeled effects. In the 

ABLRS X = 6x^ + 5x,, + fix,, . 

G V W 

Let us suppose that we have measurements {Zj^} , k-1 , ••• , 

Z|^ = H(k)x(k) + B(k)6p + V|^ , 

and an estimate x(tjt) of x(t) based on x(OjO) . 

Then 
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« X(jl|5l) + Sp(X,li)6p 
The propagation of x(a) is given by 


where 


x(Ul) = <j)X(i!.) + u(a) , 


Therefore, 


Eu(t) - 0 
Eu(ji)u{k) = . 


so that 


x(t+l jt) = 


and 


Sp(ul j&) = iSp(A|A) 


p(a+11ji) = Ex(ul l^i)x’^(.ft+l 1&) = <|)P(Ai &)<()' +q . 




V/ 




«. I 

(A" 


Now when a measurement update takes place, 


x(j^+11j!-+ 1) = x(t+l|t) + K(z^^i-Hx(a+1 |£)).: 


then 


6X(jl+ljA+l) = 6X(U1U) - KH6x(A+‘liA) 


+ KB6p + KV 


= (l-KH)[x(jl+llA)+Sp(jl+l|t)Sp3 
+ KB6p + Kv . 
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Using this, we can calculate that 

B(Jl+llJl+l) « (I-KH)P(Jl+lU)(I-KH)^ + KRK^ 


and 



S ()l+lli!,+l) = (I-KH)S (a+l|A) + KB . 

r r 

To obtain these results, we have assumed that sa , Vj^ , Wj^ and sp 
are uncorrelated random variabUs. 

Finally, the errors caused by the unadjusted parameters, 6P , are 
obtained from 


Pj -- Sp(T|T)E^'p«p)SpT(TlT) . 
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3.0 RESULTS 


The results presented in this section have all been obtained from 
computer runs implementing a 15 station grid. In base run the stations 
were placed in a 40 x 80 km area, with 20 km separation. In subsequent 
runs, grid size was varied, asd the results of this variation will be 
discussed later in this section. We have concentrated on ten specific 
effects and each is described separately. Tables have been included as 
an aid to understanding the results. A brief explanation of the tables 
follows. 

The first column of each table contains the station coordinates in 
the baseline system. In this system, the axes are different for each 
station. One axis is defined as the vector between the reference sta- 
tion, or master station, and a given station. In all of our examples, 
station 1 is the master station. The component along this axis is the 
"along" baseline eomponent, or the "L" component. Perpendicular to this 
axis in the ground plane is the second axis, where the cross component, 
or "C" component, is measured. Orthogonal to these axes in the vertical 
direction, is the height, or "H" component. Thus, the La notation, for 
example, means the along baseline component between station 2 and the 
master station. 

Also in the first column are the system biases associated with laser 
beams 1 through the number of beams. 

The recovery errors listed in the tables are in units of centi- 
meters.- A zero in any component of a station signifies that a coordinate 
system constraint was placed on that component. Coordinate system con- 
straints are discussed in Section 2.1.1 of this report. 

To the side of the table, other pertinent information is listed. 

The number of measurements is given where measurement is taken to mean 
the number of aircraft positions or the number of laser pulses. The 
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actual number of measurement‘s Is this number times the number of beams. 
Also given is the number of beams, the time Increment between pulses, and 
the aircraft low altitude. 

3,1 A ircraft Altitude and Flight Path 

One of the first priorities in the study is to determine optimal 
aircraft aUitude and flight path. Early studies showed that it is 
necessary to fly the aircraft at more than one altitude in order to 
achieve good station recovery, especially in the "height" component. Two 
altitudes seem to be sufficient in providing the necessary geometry for 
acceptable recovery. The high altitude has been chosen as 19 km, which 
is near the maximum ceiling of the aircraft. The aircraft flies faster 
and for longer periods of time at higher altitudes - two desirable traits 
in a time and cost limi?’?d mission. The low altitude controls retro- 
reflector separation - xhe lower the aircraft, the closer the stations. 
Accuracy is improved at low altitudes. Several computer runs have been 
made varying the low altitude of the aircraft to compute sensitivities. 

Table 1 shows the comparison between low altitudes of 11, 13 and 15 
km, with station separation of 20 km. If the results are compared 
directly, the picture is very confusing. Consider the 11 km and 13 km 
results for instance. For baselines closer than 45° to the flight path 
results are improved at the higher altitude. This is anomalous for two 
reasons; first, we anticipati that distance accuracy should be iso- 
tropic; second, we expect errors to increase as altitude separation 
decreases. The primary suspect for lack of symmetry is the flight path 
arrangement, particularly the lack of any A/C positions outside the long 
boundaries of the grid. 

While the higher altitude has more measurements (13 km has 440) than 
the lower altitude (11 km has 370), the improvement seems to be somewhat 
better than could -be expected on this basis alone. We hypothesize that 
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this may be because of the Increased number of measurements outside the 
grid. 


In comparing the results of the 15 km run with both the 11 km and 13 
km results> it is obvious that accuracy is lost as the two aircraft 
passes are placed closer together. 


The aircraft flies three passes at the low altitude and three passes 
at the high altitude. It is planned that these passes should be directly 
over the line of retros. However, the low pass over the boundary line 
has some difficulty in seeing six stations at all times. For this 
reason, the low passes over the boundaries were moved inside the grid. 
Results showed, however, that good viewing geometry depended strongly 
upon construction of a viewing triangle over the grid, requiring that the 
passes lie in a plane containing the retros. Therefore, both the upper 
and lower passes were moved inside, as shown in Figure^3.1-1. Table 2 
shows the comparison between flying directly over the retro line and 
flying inside the boundaries, In general, recovery is improved when 

' I 

flying inside and the improvement is probably caused by the increased 
number of measurements. 



\ 



11-15 6-10 row of targets 1-5 


Figure 3,1-1 Crossection 
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3.2 Retro Separation 


nr 



The problem of retro separation Is closely related to A/C altitude, 

At low altitudes, the retros must be closer together. Otherwise, visl- ^ 

blllty problems occur. Costwise, fewer stations placed farther apart are i 

\ 

desirable. In Table 3, a comparison Is made between two runs different 
only In the retro separation. In both runs, the low A/C lilt'^tude was 15 
km. In the first run the retros were 27.27 km apart (the largest accept- 
able separation) and In the second, they were 20 km apart. The results 
show clearly that the Interstation baselines are more accurately deter- 
mined with the smaller spacing. However, It appears that the error In 
determining a given baseline length, e.g., 80 km. Is independent of 
station spacing. 

Examination of Tables 1 and 3 In the light of this observation Indi- 
cates that a careful evaluation of A/C low altitude between U km and 15 
km, together with maximized retro spacing may produce an optimized cost 
for a given area. This will be studied in the near future. 

3.3 Number of Beams 


Another area of Interest Is the effect of varying the number of 
beams on board the A/C and varying the minimum number of retros that must 
be visible at each time for measurement processing to occur. Table 4 
shows the results of this testing. 

In column 1, the base run, 6 beams were employed along with the 
criterion that at least 6 retros must be visible at all times. In this 
run 185 A/C positions were used. In comparing these results with column 
2, where there v/ere 7 beams/6 required (185 A/C positions), we see that 
the results were somewhat better with seven beams - but improvements were 
less than 9% in baseline and generally less at long distances. In fact, 
baselines 11, 12 and 13 deteriorated slightly, apparently a numerical 
problem. Improvement was uniform In height and cross. 


3-4 


/Jt;S7MSSS I.VO TE€IL\OL(HjiaiL SYSTEMS, LSC 


In column 3, where there were 7 beems/7 required (150 A/C positions) 
results were worse. This is attributed to the fact that very few 
measurements were taken at the low altitude because of severe visibility 
problems. 

In column 4 (216 A/C positions) with 5 beams/5 visible there was 
considerable degradation, as much as Z0% on baselines 6 and 11. While 
this distinctive change in the 1, 6, 11 boundary is probably somehow 
caused by the "choose" algorithm, it is fairly clear that the five beam 
possibility is used at considerable cost in accuracy. 

Finally, column 5 employed ' .^eams/5 visible. In this run, 216 A/C 
positions were used, and results v/ere uniformly improved in baseline, 
cross and height components. Although the improvements were not great 
(generally 5-8%), the addition of these extra measurements was made with 
very little extra computational burden. 

3.4 Choose Algorithm 

The problem arises, when more -stations are visible than there are 
beams, of which stations to choose to take measurements on. Tests have 
been performed on two methods of choosing these stations. In method 1, 
the criterion for selection is distance from the A/C to the retro. 
Presently, we choose the 2 closest and 4 farthest retros in the 6 beam 
case. Method 2 utilizes a random choose algorithm. 

Contrary to initial speculation, using the random choose algorithm 
in the base run produced slightly less accurate results in all baselines 
except 7 and 12. The reason for this is not known and is thought to be 
numerical. These results may be seen in the first 2 columns of Table 5. 

In the five beam case (columns 4, 5 of Table 5) a majority of the 
baseline recoveries were better using the random choose algorithm, 
especially the longer baselines. The cross component of the longer 
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baselines seemed to improve with the random choose process, some as much 
as 1Z%, The height component of the Interior retros improved with the 
random choose, while the height recovery of the border retros deterior- 
ated with that method. 

Additional runs to see what benefit is obtained from using random 
choose will be performed, 

The 2 closest, 4 farthest criterion was originally chosen in an 
attempt to maximize geometric strength* Two characteristics of this pro- 
cedure led to the random algorithm, Fi^iSt, the closest always went to 
beam §\, which we thought led to the poorer detenai nation of bias one. 
Inasmuch as the random choose corrects this inequality, the procedure is 
effective. However, there is no significant improvement of the baseline 
recovery. Second, in a time interval when more than six stations are 
visible, the detenninistic procedure will always select the same sta- 
tions. It seemed that taking different combinations of stations might 
strengthen the solution. 

It appears, however, that any algorithm which takes both near and 
far stations will produce quite similar results. 

In an effort to support the theory of needing both close and far 
stations, a computer run was made where the 6 closest stations were 
always chosen. The results of this run can be seen in column 3 of Table 
5. Clearly, the baseline recoveries have suffered as a result of poor 
geometry. 

3 . 5 Basic Mea surement Noise Variance Versus Number of Measurements 


Computational burden is greatly increased by a decrease in time 
between measurements. When the At is halved, central processing unit 
time is doubled. For this reason, we were interested in performing tests 
to find out if a large number of measurements could be simulated by using 
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a smaller basic measurement noise variance. The results of this test can 
be seen in Table 6 where the baseline errors are virtually identical. In 
column 1 these errors were achieved using a At » 2.5 and a basic 
measurement noise variance = 1. In column 2, At «s 5.0 and the basic 
measurement noise variance = 0.5. 

It is interesting to note in Table 7, columns 1, 2, 3 (At = 10, 5, 
2.5) that behavior of the cross component is very well described by 

. This same behavior is shown in columns 4 through 8 (measurement 
o « 1.414, 1.0, 0.707, 0.5, 0.3). We see also that halving At is 
equivalent to halving measurement noise variance. 

This indicates that the observability problem, or a priori influ- 
ence, which strongly affects baseline and height is not impacting the 
cross component. 

3.6 A Priori Variance 

In Table 8 we decreased the a priori variance on the bias. In 
column. 1 the variance was 1, in column 2 the variance was 0.25 and in 
column 3 the variance was 0.09. These runs indicate that prior informa- 
tion on bias is still strongly affecting height and baseline recovery. 

The cross component, however, is independent of bias. Previous runs have 
shown that recovery is independent of prior data on station and A/C - at 
least while that prior is at the 10m - 30m level. 

These results are consistent with those in Table 7 indicating that 
the cross component is noise limited and independent of bias. In order 
to determine whether the difficulties in reaching the noise-limited con- 
dition are because of marginal observabil ity or the influence of prior 
data, two runs were made with large prior variances on bias. These runs, 
column 4 (obias = 1000 cm) and column 5 (ubias = 500 an), show that 
results are Independent of prior variances when the biases are much 
larger than the posterior estimates. This indicates that the system is 
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completely observable. However, It does not tell us when the system 
becomes noise limited. In order to detennine how large numbers of 
measurements will affect the recovery, the simplest procedure is to 
decrease the measurement noise and properly adjust the process noise. 

3.7 Effect of Process Noise 


The results presented in the previous sections have all assumed A/C 
position uncertainty of 9.5 m in each axis at each measurement point. 
Since the position of the A/C must be solved for and the uncertainty 
cannot grow to this magnitude in the 0.1 - 0.2 second interval between 
"shots", this uncertainty induces a very pessimistic estimate of recovery 
errors. To provide a realistic evaluation, a process noise simulation 
has been implemented and is described in Appendix A. In order to check 
She implementation, unreasonably large values were used. Tliese results 
are shown in Table 9. 

An estimate of how large the uncertainties should be is difficult to 
obtain. However, using the parameters described in Appendix A and a 10 
sec measurement interval, there is no improvement. Hov-zever, if the 
interval were reduced, as it is in the operational system, improvement 
would show. It is difficult to evaluate this properly without actually 
running at 0.2 sec. however, which is financially prohibitive. 

In the operational data reduction, the process noise can be evalu- 
ated, of course, so as to close the loop on the numerical values. 

3.8 Refraction 


Since refraction is a very important error source, it is necessary 
to consider its contribution to the errors in baseline recovery. Table 
11, column 2 shows our first attempt at refraction in the "consider" 
mode. The noise only errors are identical to the base run errors in 
column 1. The unadjusted errors are all ,ess than 0.7 cm. Column 3 
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shows the results of the same run when the random choose algorithm is 
selected. The noise only baseline errors exhibit the same character- 
istics as previously discussed, while the unadjusted errors show "random" 
improvement throughout the grid. Column 5 shows the results when the 
time increment between measurements was reduced to 5 sec. Hie noise only 
errors are identical to the errors in the base run for At « 5 sec. 
(column 4). 

3.9 Cross Passes 

Four cross passes were added to the A/C flight path to see if the 
geometry would be improved enough for significant reductions in the base- 
line recoveries. The results in Table 12 . show that the baseline com- 
ponent errors were slightly better than would have been expected given 
the additional number of measurements and that the height component 
errors were very close to the expected errors. 

Therefore, it seems that there is no pronounced improvement in 
recoveries as a result of adding cross passes to the flight plan. 

3.10 Operational Projections 

From Table 11, an accurate estimate of refraction errors is 
obtained. From Table 7, the evolution of error can be traced as At 
becomes smaller, and it appears that each halving of At is accompanied 
by an error reduction of about IBJS. This reduction is applied six times 
to obtain the equivalent of a pulse interval of 0.156 sec. Although the 
laser pulses every 0.1 sec we anticipate a 33fi signal loss and thus 0.15 
sec seems a reasonable data interval. These results are multiplied by 
(to simulate the errors in baseline change obtained from two experi- 
ments) and are shown in Figure 3.10-1, 

In Figure 3.10-2 the analogous baseline rate accuracy is shown. 
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While larger grids have not been tested, these results v/ill hold for 
the same size subsections of larger grids, provided the number of 
measurements and the geometry is preserved. 





























Tabla 4 - Variation In Numbar of Beams 


Stations 


B I asos 



6 beams 
6 visible 

7 beams 
6 visible 

7 beams 
7 visible 

5 beams 
5 visible 

6 beams 
5 visible 

H2 

1.51 

1.41 

1.68 

1. 58 

1.39 

L2 

0.77 

0.71 

0.83 

0,85 

0.72 

C2 

a 51 

0.46 

0,51 

a?i 

a50 

H3 

1.69 

1.77 

2. 16 

1.95 

1.73 

L3 

1. 13 

1. 04 

1,26 

1. 19 

1,04 

C3 

a 51 

0.49 

0.52 

0.70 

a 50 

1*4 

1.53 

1.42 

1.69 

1.58 

1.40 

U4 

1.55 

1,42 

1.73 

1. 63 

1.42 

04 

0.55 

a46 

0.52 

0.81 

a 54 

HS 

0.00 

a 00 

0.00 

0,00 

0.00 

LS 

2.05 

1.87 

2.28 

2. 12 

1,86 

C5 

0.00 

a 00 

0.00 

a 00 

0.00 

H6 

1.62 

1. 59 

1.86 

1.81 

1. 50 

14 

0.82 

0.79 

0,92 

1.06 

a79 

C6 

0.73 

a 70 

a8i 

0.86 

a 70 

H? 

1. 07 

1.12 

1. 18 

1.26 

1,05 

L? 

0.95 

a9i 

0.97 

1.09 

a 92 

c? 

0.75 

1 a 68 

0.92 

0.89 

a 70 

Ha 

I 1. 15 

1. 15 

1.30 

1.35 

1. 10 

Lb 

1.20 

U16 

1.33 

1,33 

1. 14 

ca 

0.79 

0.74 

a 97 

a 97 

0.74 

H9 

1. 11 

1.01 

1, 18 

1.41 

1.07 

L9 

1.54 

1.46 

1. 73 

1.68 ' 

1.45 

C9 

0.63 

a 78 

a 98 

1.05 

0.78 ■ 

Hio 

1.60 

1.48 

1.83 

1.92 

1. 53 

Lio 

1.98 

1.85 

2.23 

2.12 

1.84 

CIO 

a 87 

a 87 

1. 02 

1. 19 

0.63 

Hll 

2,86 

2.65 

3.32 

3.27 

2.65 

Ln 

1.23 

1. 24 

1.38 

1.65 

1.21 

Cii 

1. 20 

1. 14 

1.55 

1.41 

I. 15 

Hi 2 

1.24 

1. 16 

1.43 

1.46 

1. 16 

L12 

1.28 

1.30 

1.37 

1. 63 

1.26 

Cl 2 

1. 13 

1.04 

1.36 

1.43 

1.09 

H13 

0.00 

a 00 

0.00 

0,00 

0.00 

L13 

1.43 

1.44 

1.54 

1.70 

1.40 

Ci3 

1.21 

1. 10 

1.46 

1.54 

1. 14 

Hi 4 

1,21 

1.19 

1.46 

1.42 

1. 15 

Ll4 

1. 70 

1.67 

1.88 

1.96 

1.63 

Cl 4 

1.25 

1. 17 

1.54 

1.64 

1.19 

Hi 5 

2.76 

2.70 

3.28 

3. 16 

2,62 

LlS 

2,08 

2.01 

2.33 

2.33 

1. 97 

CIS 

1.32 

1.26 

1.66 

1.74 

1.24 

1 

0.43 

0,41 

a 43 

0.49 

a 42 


0,40 

a 38 

0.39 

a44 

a 40 

3 

0.42 

a39 

a4i 

0.48 

0.41 

4 

0,41 

a 38 

a39 

a45 

a 40 

5 

a40 

a 37 

0.38 

0.45 

0.39 

6 

0.39 

0,37 

a 38 

- 

0.39 

7 

- 

a 36 

0.37 


- 


Low Al titudo ■ 11 km, 
At » 10. soc. 


Co I umn 


Number of 
Moasuremonts 
Per Beam 


1 165 


2 165 

3 150 


4 216 


5 


216 











Table 6 - Belle Heosuresent Noise Vorlence 
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Tobla 7 


Stations H2 
LZ 
02 


At»IO. Ut»5. At*-2.5 


0,82 0.63 

1.26 1.07 

0.59 0.41 


2.11 1.57 

0. 94 0. 74 

0.86 0.60 


m»as. <t moos. 0 moos. 0 moos, 0 moos. 0 
*1.414 ■ 1, " .707 ".5 " ,3 


Low Altitude ■ 11 km 

6 Booms 

At m 5 , soc. 

<Col. 4-8) 


Column 

Number of 
i4oo,suroffl«nts 
Par Beam 

1 

185 

2 

370 

3 

738 

4 

370 

5 

370 

6 

370 

7 

370 

s 

370 


Biases 



0.40 

0.37 

0.39 

0.39 

0.38 

0,37 












Table 9 - Process llolso Chockout 
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Tablo 10 - Process Noise 


Stations 


Biases 



Process Noise 
Not Model led 

Process Noise 
Mode 1 led 

H2 

1.51 

1. 51 

L2 

0.77 

0,77 

C2 

0,51 

0. 51 

H3 

1.89 

1.90 

1-3 

1. 13 

1, 13 

C3 

0, 51 

0.51 


1.53 

1. 53 


1.55 

1.55 

O 4 

0.55 

0.55 

Hs 

0.00 

0.00 

L5 

2.05 

2.05 

C5 

0.00 

0.00 

He 

1.62 

1.62 

Le 

0.82 

0.82 

C6 

0.73 

0.73 

H? 

1.07 

1.07 

L? 

0.95 

0,95 

C7 

0.75 

0.75 

Ha 

1.15 

1. 15 

La 

1.20 

1.20 

c« 

0.79 

0.79 

H9 

1. 11 

1.11 

L9 

1.54 

1.54 

C9 

0.83 

0.83 

HlO 

1.60 

1.60 

LlO 

1.98 

1.98 

CIO 

0.87 

0.87 

Hll 

2.88 

2. 88 

Lll 

1.23 

1.23 

Cll 

1.20 

1.20 

Hi 2 

1.24 

1.24 

L12 

1.28 

1.28 

Cl 2 

1. 13 

1. 14 

H13 

0.00 

0.00 

L13 

1.43 

1.43 

Cl 3 

1.21 

1.21 

Hm 

1.21 

1.21 

Lm 

1.70 

1.70 

Cl4 

1.25 

1.25 

His 

2.76 

2.76 

LlS 

2.08 

2.09 

CIS 

1.32 

1.32 

1 

0.43 

0.43 

2 

0.40 

0.40 

3 

0.42 

0.42 

4 

0.41 

0.41 

3 

0.40 

0.40 

6 

0.39 

0.39 

— 


Low Altitude ■ 11 kffl, 
6 Beams 


At >• 10 

sec. 


Number of 

Column 

Measurements 


Per Beam 

1 

185 

2 

185 


3-21 
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Toblo 11 - Ro fraction 


Stations H2 

C2 


Ho I so Uncd- Ho I so Unad- 
Only Juetod Only Justod 

1.51 0.28* 1.62 0.16 

0.77 0. 14 0.79 0. 12 

0.51 0.03 0.54 0.04 


4» 5» 

Nolso Unad- 
Only Justod 

1.16 l.''l6 0.31 

0. 59 0. 59 0. 1 5 

0,36 0,36 0.02 

1.47 1.47 0.43 

0.91 0,91 0.32 

0.36 0.36 0,02 

1.17 1.17 0.34 

1.20 1.28 0.49 

0,39 0.39 0.02 


Biases 



Low Altitude » 1 1 km. 

6 Beams 

Cols, 1-3 At » 10 soc. 
Cols, 4-5 At » 5 soc. 


Column 

Number of 
Measi/rements 
Po'A Beam 

1 

185 

2 

185 

3 

185 

4 

370 

5 

370 


I» Base Run - 
Refraction Nbt 
Considered 

2* Base Run - 
Refraction 
Considered 

34 Random Choose - 
Ro fraction 
Considered 

4*^ Base Run - 
Refraction Hot 
Considered 

5* Base Run - 
Refraction 
Considered 





Table 12 - Effect of Cross Posses 
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4.0 Conclusions 


Figure 3.10-1 shows clearly that centimeter accuracy is achieved 
over 80 km baselines, even without exotic refraction equipment - such as 
two-color lasers. 

It appears (Section 3.3) that the laser must be provided with six 
independent beams. Current investigations assume these to be indi- 
vidually steerable. 

To minimize retroreflector cost while obtaining maximum grid extent, 
there appears to be an optimum aircraft low altitude between 11 km and 15 
km and an associated retro spacing between 20 km and 27 km. 
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Appendix A OF POOR QUALITY 

COMPUTATION OF A/C STATISTICS 

The aircraft is assumed to be moving in a straight line during each 
pass, this is the reference trajectory. However, there are perturbations 
in the flight path which cannot be obtained from the onboard navigation 
or from GPS; these constitute the uncertainties in the aircraft position 
and velocity. We assume that these perturbations are generated by a 
second order system forced by white noise. Because these perturbations 
are deviations from the reference path, we assume that there is some con- 
trol system which drives the aircraft back to the reference path. This 
provides a rationale for making the linear system asymptotically stable 
and thereby making the perturbation sequence stationary. 

The model we have chosen is 


I p + ZXp + (X-+v^) p = u 

" with 

Eu = 0 , Eu(t) u(t) = q«S(t-T) . 

I- 

This model is physically inaccurate, since the perturbing accelerations, 
u , will not be white. However, we expect this to portray a satisfactory 
representation, on the basis that the combination of pilot/aircraft 
dynamics and atmospheric turbulence will generate a position sequence 
having this kind of correlation. 

f 

i 

’ Initially, we shall presume that there are three independent Markov 

processes (x,y,z) with identical characteristics, while this may be 
unrealistic, it will be more realistic than the current procedure and 
will enable us to evaluate the benefits of this approach. In any case 

I 

I the program can be modified easily to provide for different character- 

i sties. 
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The program currently propagates the A/C covariance via the usual 

P(t+t) ^ 4.p(t) i ^ (K, (1) 

when A / 0 . When ^ » 0 , however, 

P't">A/C “ PoA/C 

where P^^^^ is the initial aircraft covariance. In order to change 
this, we assume that x, y, and z each obeys a system 



0 1 ^ 
-(X^+v^) - 2X 



u , 


I 


so that 


<Ji r e 


-xt 


cos vt + — sin vt 


X^+v^ 


sin 


sin vt 


cos vt - -- sin vt 


and 


where 


'^11 ^12’ 

^12 ^22 


^11 = 


-2Xt 


4X(X^+v^) 4v' 


1 X cos 2vt - V sin 2vt 
■ X"+v 2 


•-li: 


-2Xt 

5 — (1 - cos 2vt) 

4v^ 
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o » 1 X e 

^22 V: 


-2Xt 


[- 


+ ~ cos + s1n 2 vt 


]■ 


and q is the spectral density of u . 
In equilibrium (t = , we have 


4X(x^+v^) 

qj2 “ 0 

^22 “k- 


These values are derived in the appendix. 

The parameter values have been selected to approximate A/C dynamic 
characteristics with 

10 sec 
and 

V = - ■ 2^ 

30 sec ’ 


which gives 


qjj M = 46.41 
qj2 (") = 

422 <“) = 2.5 . 

We want the equilibrium value to have = 9E6 . Therefore 
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q » 194000 ■2— and * 485000 - — r* (696 cm/sec)^ , a very large 
sec^ ^ sec 

uncertainty. 

The way this Is handled 1n the computer, Is that Q Is computed for 
each step and the A/G uncertainty updated according to equation (l). 
Whenever a new track is started, however, the A/C uncertainty returns to 
. To complicate matters* we have noticed In the past that when 
station uncertainties are small and A/C uncertainties are large then the 
matrix 


CHPh’^ + R] , 

appearing In the covariance update, Is numerically singular. This occurs 
when the uncertainties in A/C position are so much larger than noise or 
station uncertainties that the matrix appears to have rank less than the 
number of beams. This difficulty has been handled by processing the 
measurements one at a time. 
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Appendix B 

Covariance Matrix of a Second Order Markov Process 


Consider the system 


^ “ [^-(XW) -2xJ ’ 

< + 

[:]■ 

( 1 ) 

associated with the equation 




X + 2X X + (X^+v^) ; 

< “ u . 

(2) 

The solution of this system is 




x(t) = 4'(t)x(0) + / 4>(t-s) 
0 

~Q~ 

1 

u(s)ds , 

( 3 ) 

x(t) = *(t) x (0) + w . 

► 

( 4 ) 

If we take u to be a Gaussian random process, then it is 
reasonable to discuss the covariance of x . Let 


E u(t) = 0 , E u(t)u(T) 

= q 

6(t-t) 

( 5 ) 


where q is the (constant) spectral density of u , and define 

P(t) = E x(t) x'^(t) . (6) 

Then 

P(ttx) = H^) P(t) »^(T) + Q(t) q . (7) 
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Expl Icitly, 


^(t) » 


cos Xt - sin vt 




sin vt 


sin vt 


cos t - sin vt 


We now need to compute Q(t) , 


Q(t) « E w(t) w'(t)/q . 


This can be obtained directly from 


w(t) « / 


sin v(t-s) 

V 

cos v(t-s) - ^ sin v(t-s) 


( 8 ) 


(9) 


u(s)ds ; (10) 


which comes from equation (3). Using (1) and (5) it is clear that 
E w(t) = 0; and 


qii 


(t) = 1- / sin^ v(t-s) ds , 


0 


( 11 ) 


912 ^ ^ sin v(t-s) cos v(t-s) ds - Xq (t) , 

0 


( 12 ) 


and 


(■^) = / cos^ v(-r-s) 


ds - 2Xq^^(^) - ^Sii(T) . (13) 


These can be rewritten as 


qn (’') = •^ / e 


-?Xs 1 - cos 2VS j_ 
2 ds . 
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^12 



] -2Xs sin 2\»s 
J e ■ p 


ds " * 


and 


^22 (0 = 1 ds - 2Xq^2(T) - X^„(T) 

from which the values can be calculated: 

1 n X cos 2^>T .. V sin 2 vt 


''i:'’'’ ' ;x(xW) 


“1 \ COS 2^>T - V sin 

?1T~ T “ ^^2 ^ ^2 J 


qiaCO 


I X sin 2v-t + ^^cos 2v.l _ Xq^^(T) 


4(X^+v^) 4v 


-2 Xt 

® (1 - cos 2\>t) , 


X^ + 


4v' 


”22 <") “ 


2X^+v^ 

4X(X^+v^) 


-2Xt 


'1 ^ X cos 2VT - V sin 2 vt 


X^ + 


] 


2Xq^2 ^11 


, -2 Xt 

1 + e - 

+ — J- 


[ 


x^+v^ i L- 
Xv^ 


cos 2 vt 


; sin 2^'^'^ 


These formulae have been checked by verifying that 

q^(0) = 0 = q^2^°^ ^22^°^ 


and 


= 0 = qi2^°^ ’ ^22^^^ ” ^ 

which agrees with equations (U) - (13). 
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